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Infrared-terahertz  double-resonance  spectroscopy  of  CH3F  and  CH3C1  at  atmospheric  pressure 
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A  method  for  highly  selective  remote  sensing  of  atmospheric  trace  polar  molecular  gases  is  described.  Based 
on  infrared-terahertz  double-resonance  spectroscopic  techniques,  the  molecule-specific  coincidence  between 
the  lines  of  a  CO2  laser  and  rotational-vibrational  molecular  absorption  transitions  provide  two  dimensions 
of  recognition  specificity:  infrared  coincidence  frequency  and  the  corresponding  terahertz  frequency  whose 
absorption  strength  is  modulated  by  the  laser.  Atmospheric  pressure  broadening  expands  the  molecular  recognition 
“specificity  matrix”  by  simultaneously  relaxing  the  infrared  coincidence  requirement  and  strengthening  the 
corresponding  terahertz  signature.  Representative  double-resonance  spectra  are  calculated  for  prototypical 
molecules  CH3F  and  CH3C1  and  their  principal  isotopomers  from  which  a  heuristic  model  is  developed  to 
estimate  the  specificity  matrix  and  double-resonance  signature  strength  for  any  polar  molecule. 

DOI:  10. 1 103/PhysRevA.85. 052507  PACS  number(s):  33.20.-t,  33.40.+f 


I.  INTRODUCTION 

Double-resonance  (DR)  techniques  are  widely  used  in  all 
phases  of  optical  spectroscopy,  from  ultrafast  pump-probe 
and  excitation  correlation  measurements  in  the  ultraviolet, 
visible,  and  infrared  (IR)  spectral  regions  to  investigations  of 
molecular  collisional  physics  in  the  microwave-microwave  and 
infrared-millimeter  wave  regions  [1-7],  Although  wavelength- 
degenerate  DR  techniques  provide  some  insight  into  the  relax¬ 
ation  dynamics  from  a  photoexcited  state,  wavelength-non- 
degenerate  DR  techniques — which  rely  on  the  ability  to  tune 
the  frequencies  of  the  pump  and  probe  beams  independently — 
provide  greater  insight  into  energy-transfer  pathways  and  rates 
[8-11].  Sometimes,  the  limited  frequency  tunability  of  pump 
or  probe  sources  limits  the  effectiveness  of  DR  techniques.  For 
example,  optically  pumped  far-infrared  lasers — which  use  a 
line-tunable  CO2  IR  laser  to  excite  a  terahertz  (THz)  rotational 
population  inversion  in  a  low-pressure  gas — are  constrained 
by  the  rare  coincidence  (<100  MHz)  between  any  CO2  laser 
line  and  any  molecular  rotational-vibrational  IR  transition 
[12-14].  This  rare  coincidence  forces  optically  pumped  far-IR 
(OPF1R)  lasers  to  operate  on  only  a  few  THz  frequencies 
per  molecule,  but  as  the  gas  pressure  and  corresponding  IR 
and  THz  linewidths  increase,  greater  tunability  and  spectral 
coverage  are  possible  [15-20]. 

Recently,  we  have  proposed  that  the  IR-THz  DR  technique 
may  be  adapted  for  remote  sensing  of  trace  gases  in  the 
atmosphere  if  IR  pulses  of  sufficiently  short  duration  (~  100  ps) 
can  be  generated  [21].  In  particular,  it  was  shown  that  a 
100-m-thick  cloud  of  CH3F  with  a  uniform  concentration 
of  1  ppm  (~  1  -mTorr  partial  pressure)  could  be  detected 
using  a  THz  transceiver  tuned  to  the  frequency  of  the  IR 
pump-induced  THz  DR  signature.  Detection  occurs  through 
the  co-propagation  of  the  IR  pump  and  the  THz  probe  beams 
through  the  trace  gas  cloud  to  a  retroreflecting  surface  up 
to  1  km  away.  If  absorbed  by  the  trace  gas,  the  pulsed 
IR  pump  beam  briefly  alters  the  absorption  strength  of  the 
returned  continuous-wave  THz  probe  beam,  and  this  repeating 
modulation  may  be  detected  at  the  pulse  repetition  frequency 


of  the  laser.  The  actual  sensitivity  of  this  technique  depends  in 
a  complex  manner  upon  the  IR  and  THz  absorption  coefficients 
of  the  molecule  and  the  atmospheric  water  vapor  content. 
Recognition  specificity  is  achieved  through  the  molecule- 
unique  combination  of  IR  pump  coincidence  frequency  and  the 
frequency  of  the  pump-induced  change  in  the  THz  signature. 
For  a  given  molecule,  only  a  few  of  its  hundreds  of  IR 
rovibrational  transitions  are  coincident  with  any  CO2  laser 
line,  and  each  rare  coincidence  produces  a  unique  pump- 
induced  absorption  change  in  a  few  of  the  hundreds  of  THz 
rotational  transition  frequencies.  These  two  dimensions  of 
recognition  (IR  pump  and  THz  DR  frequencies)  comprise  a 
“specificity  matrix”  whose  sparseness  can  be  used  to  identify 
an  atmospheric  trace  gas  remotely  with  enough  spectral 
discrimination  that  even  isotopomers  may  be  distinguished. 

One  of  the  findings  of  that  paper  [21]  is  that  the  coincidence 
requirement  between  CO2  laser  and  IR  molecular  transition  is 
relaxed  when  the  trace  gas  primarily  collides  with  atmospheric 
pressure  gases  (see  Fig.  1).  Specifically,  Fig.  2  indicates 
that  the  rovibrational  IR  linewidth  of  CH3F  evolves  from 
a  Doppler-broadened  Avp  =  33  MHz  [half  width  at  half 
maximum  (HWHM)]  for  low-pressure  operation  (typical  of 
OPFIR  lasers)  to  a  pressure  broadened  A Vp  =  2.3-GHz 
HWHM  when  the  same  1  mTorr  of  gas  is  in  an  atmosphere  of 
N2  and  O2  [22],  Moreover,  the  IR  linewidth  function  changes 
from  Gaussian  to  Lorentzian,  meaning  that  the  spectral  overlap 
falls  off  less  abruptly  with  frequency  at  atmospheric  pressure. 
This  broadening  relaxes  the  pump  coincidence  requirement 
from  those  rovibrational  lines  within  66  MHz  of  a  CO2  laser 
line  [spaced  24-69  GHz  apart,  see  Fig.  1(a)]  to  those  within  4.6 
GHz  [23].  The  benefit  is  that  favorable  IR-THz  pump-probe 
combinations,  including  transitions  that  are  not  accessible 
at  low  pressure,  may  be  selected  for  atmospheric  remote 
sensing. 

In  order  to  explore  the  molecule-unique  signatures  that 
constitute  the  recognition  specificity  matrix,  we  calculate 
the  double-resonance  signatures  of  all  IR-THz  pump-probe 
coincidences  at  atmospheric  pressure  for  12CH3F  and  CH3  35C1 
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FIG.  1.  Infrared  spectra  of  (a)  C03  laser  [23],  (b)  CH3F,  and 
(c)  CH3C1  [28,38]  between  900  and  1100  cm'1.  Rotational  THz 
spectra  of  (d)  J  =  8  to  9  in  V3  =  1  of  CH3F  and  (e)  J  =  8  to  9  in 
V6=  1  of  CH3C1. 


and  their  principal  isotopomers  13CH3F  and  CH3  37C1  that 
are  1%  and  24%  abundant,  respectively.  These  molecules 
were  chosen  as  prototypical  because  their  excited  vibrational 
levels  are  well  isolated  but  different  in  nature:  A  CCU  laser 
that  is  line  tunable  over  a  9-1  l-/tm  wavelength  region  (900- 
1100  cm-1)  excites  the  V3  stretching  mode  of  CH3F  and  the  V( 
bending  mode  of  CH3CI.  We  show  that  the  DR  signatures  are 
dramatically  different  for  these  two  types  of  vibrational  states 
in  ways  that  are  unique  to  the  molecule,  providing  additional 


FIG.  2.  (Color  online)  The  calculated  n(v^)  for  each  <2a:(12) 
transition  of  12CH3F  pumped  by  the  9P(20)  CO2  laser  line,  overlaid 
by  the  HWHM  linewidth  of  IR  rovibrational  transitions  of  continuous 
lines:  CH3F,  and  dashed  lines:  CH3C1,  from  1  to  760  Torr. 


recognition  specificity.  Because  many  new  coincidences  occur 
at  atmospheric  pressure,  we  show  how  to  estimate  the  strength 
of  the  DR  signals  so  that  the  IR-THz  combinations  most  useful 
for  remote  sensing  may  be  identified. 


II.  DOUBLE-RESONANCE  SPECTROSCOPY 
A.  Pump  coincidences 

The  energy  levels  of  the  ground  and  excited  vibrational 
states  of  CH3F  and  CH3CI  were  calculated  using  the  rotational 
constants  listed  in  Tables  I  and  II.  For  CH3F,  the  standard  P- 
type  (A/  =  —  1),  0-type  (A/  =  0),  and  R-type  (AJ  =  +  1) 
transitions  from  ground  state  to  the  V3  vibrational  level 
produce  a  classic  stretching  mode  spectrum  [Fig.  1(b)]  whose 
absorption  strengths  are  determined  by  a  combination  of  the 
ground  state’s  fraction  of  the  temperature-dependent  parti¬ 
tion  function,  quantum-mechanical  transition  matrix  element 
(Table  III),  and  dipole  derivative  for  the  vibrational  transition. 
The  notation  for  a  rovibrational  transition  is  X K(J ),  where 
J,K  designate  the  rotational  state  of  the  origin  in  the  ground 
vibrational  level  and  X  =  P ,  Q,  or  R  depending  on  the  type 
of  rovibrational  transition  excited. 

By  comparison,  the  spectrum  of  the  Vf,  bending  mode  is 
characterized  by  /  doubling  that  couples  the  bending  motion 
/  =  ±  1  with  the  rotation  designated  by  the  K  quantum  number. 
The  notation  is  modified  to 3  Xr  (J)  to  account  for  whether  the 
rovibrational  transition  excited  in  Vg  satisfies  Z  =  +1(Y  =  R) 
or  —1  (Y  =  P).  The  corresponding  spectrum  [Fig.  1(c)]  is 
remarkably  different  than  for  the  bending  mode  in  CH3F 
because  /  doubling  breaks  the  ±K  degeneracy,  dramatically 
affecting  the  energy-level  structure  and  transition  matrix 
elements  (Table  III).  An  examination  of  a  partial  low-pressure 
rotational  spectrum  of  these  two  excited  states  [Figs.  1(d)  and 
1(e)]  shows  the  effects  of  /  doubling  and  hyperfine  splitting 
in  CH3CI.  (Hyperfine  splittings  in  CH3CI  that  are  resolved  at 
low  pressure  may  be  ignored  at  the  atmospheric  pressures  of 
interest  here.) 

Using  these  calculated  energy  levels,  the  allowed  IR  rovi¬ 
brational  transition  frequencies  were  calculated  and  were  com¬ 
pared  with  the  known  frequencies  produced  by  line-tunable 
CO2  lasers  [23].  To  find  nearly  overlapping  coincidences, 
assume  for  the  moment  that  the  IR  pump  is  a  continuous-wave 
CO2  laser,  so  the  coincidence  occurs  within  the  linewidth  of  the 
molecular  rovibrational  transition.  For  a  pure  1-mTorr  gas  of 
CH3F  or  CH3CI  at  300  K,  the  salient  IR  rovibrational  transition 
is  Doppler  broadened  with  a  Gaussian  HWHM  linewidth  AvD 
of  33.4  and  27.4  MHz,  respectively  [24].  All  rovibrational  IR 
frequencies  are  then  calculated,  and  the  “offset”  frequencies 
V/R  —  npump  of  those  nearly  coincident  with  specific  CO2  laser 
lines  are  compared  with  the  values  known  from  extensive 
research  in  OPFIR  lasers  [25,26].  The  results,  shown  in 
Table  IV,  indicate  the  offset  frequencies  are  in  good  agreement 
when  the  vibrational  energies  in  Table  I  are  used.  As  a  notation 
to  simplify  the  discussion,  the  rotational  states  connected 
by  the  pump  are  labeled  (Jp,Kp,V 0)  and  (J'p,K'p,Vi).  The 
principal  rotational  transitions  exhibiting  THz  DR  signatures, 
involving  the  two  pairs  of  adjacent  rotational  states  ( J p  ±  1, 
Kp,V 0)  and  (Jp  ±  1  ,Kp,Vi),  are  labeled  R q  and  Rjj",  R~,  and 
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TABLE  I.  Molecular  constants  of  methyl  fluoride. 


12 

CH,F  [34] 

13CH3F  [35,36] 

0 

11 

;£■ 

V,  =  1 

0 

II 

V,  =  1 

A  (MHz) 

15  5352.70  (36)“ 

155  058.474(19)“ 

155  365.345  (40)“ 

155  076.4480  (309)“ 

B  (MHz) 

25  536.149  80  (54)“ 

25  197.510  60(87)“ 

24  862.665  (2 )“ 

24  542.131  (3)“ 

Dj  (kHz) 

60.2214(20)“ 

56.8868  (33)“ 

57.724  35  (1130)“ 

55.069  39  (1050)“ 

D]t  (kHz) 

439.8156  (240)“ 

518.2386  (471)“ 

424.83  (11)“ 

477.773  (140)“ 

Dk  (kHz) 

2108.0  (75)“ 

2014.9000  (3358)“ 

2137.293  (1881)“ 

2071.5  (39)“ 

H,  (Hz) 

-0.03214  (93)“ 

-0.193  20  (29)“ 

-0.021  30  (728)“ 

-0.1200(66)“ 

Hjjk  (Hz) 

1.9470  (275)“ 

16.076  90  (5156)“ 

1.5360  (872)“ 

9.7025  (1510)“ 

Hjkk  (Hz) 

22.3610(1601)“ 

-96.449  47  (48566)“ 

21.500  (764)“ 

-38.421  (1130)“ 

Hk  (Hz) 

0 

0.119  7871  (22754)“ 

150.4  (2471)“ 

231.99  (138)“ 

Lj  (mHz) 

0 

0.010  20  (87)“ 

0 

0 

L  jjjk  (mHz) 

0 

-1.7468  (187)“ 

0 

-0.7732  (594)“ 

L jjkk  (mHz) 

0 

29.5648  (2138)“ 

0 

13.300  (504)“ 

Ljkkk  (mHz) 

0 

-173.2365  (14660)“ 

0 

-89.64  (270)“ 

Lk  (mHz) 

0 

0 

0 

0 

Vibrational  energy  (cm-1) 

- 

1048.610  701  (10)“ 

- 

1027.49320  (2)“ 

Dipole  (debye)  [31] 

1.857 

1.857 

1.857 

1.857 

Dipole  derivative  (debye)  [31] 

- 

-0.271 

- 

-0.271 

“Uncertainties  in  the  least  significant  digit(s)  are  in  parentheses. 


R~f  ,  respectively  (see  Fig.  3).  The  frequencies  of  these  four 
rotational  transitions  are  presented  in  Table  IV. 

Next,  the  rovibrational  transitions  of  CH3F  or  CH3CI  with 
1  mTorr  of  partial  pressure  are  broadened  by  an  atmosphere  of 
N2  and  O2  (i.e.,  concentration  0  ~  1  ppm  of  analyte),  and 
the  coincidence  calculation  is  repeated.  Using  a  pressure¬ 


broadening  parameter  of  3-MHz/Torr  FWHM — typical  for 
hard-shell  collisions  with  nonpolar  molecules  [22] — the  esti¬ 
mated  linewidth  of  Au;)  =  2.3-GHz  HWHM  far  exceeds  the 
Doppler-broadened  linewidth.  Using  a  Lorentzian  line-shape 
function  with  this  linewidth  parameter,  the  overlap  expands, 
and  many  new  coincidences  are  found. 


TABLE  II.  Molecular  constants  of  methyl  chloride. 


12ch3 

35C1  [37] 

12ch3 

37C1  [37] 

0 

II 

s? 

v6  =  i 

0 

II 

v6  =  1 

A  (MHz) 

156  051.1 

156  810.71  (2 )“ 

156  440 

157  197.02  (3)“ 

B  (MHz) 

13  292.8752  (63)“ 

13  243.8840  (69)“ 

13  088.1660(93)“ 

13  039.967  (10)“ 

Dj  (kHz) 

18.1010(39)“ 

18.1360  (45)“ 

17.5630  (78)“ 

17.5960  (87)“ 

Djk  (kHz) 

198.77  (13)“ 

203.37  (14)“ 

193.48  (13)“ 

197.85  (28)“ 

Dk  (kHz) 

2653.1 

2735.2  (1)“ 

2501.0 

2583.0(1)“ 

Hj  (Hz) 

-0.008  45  (840)“ 

-0.010  10  (96)“ 

-0.0121  (20)“ 

-0.0140  (23)“ 

Hjjk  (Hz) 

0.312(33)“ 

0.375  (36)“ 

0.312(33)“ 

0 

Hjkk  (Hz) 

9.53  (66)“ 

5.16(72)“ 

8.39  (66)“ 

4.08  (190)“ 

Hk  (Hz) 

0 

0 

0 

0 

Af  (MHz) 

0 

39  275.246  (30)“ 

0 

39  560.346  (33)“ 

q  (MHz) 

0 

3.639  735  (250)“ 

0 

3.580  87  (33)“ 

q,  (kHz) 

0 

-0.017  24  (10)“ 

0 

-0.017  850(17)“ 

9jj  (Hz) 

0 

0 

0 

0 

rjj  (MHz) 

0 

0.467  059  (78)“ 

0 

0.457  777  (60)“ 

m  (MHz) 

0 

4.027  59  (100)“ 

0 

3.4039  (15)“ 

qjk  (kHz) 

0 

0.128  60  (57)“ 

0 

0.1267  (10)“ 

*ljj  (Hz) 

0 

-0.6601  (250)“ 

0 

0 

Vkk  (kHz) 

0 

-0.1020(45)“ 

0 

-0.1250  (75)“ 

Vibrational  energy  (cm-1) 

1018.071  10(7)“ 

1017.695  45  (1)“ 

Dipole  (debye)  [31] 

1.8989 

1.8989 

1.8989 

1.8989 

Dipole  derivative  (debye)  [31] 

-0.0388 

-0.0388 

“Uncertainties  in  the  least  significant  digit(s)  are  in  parentheses. 
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TABLE  III.  Branching  ratios  for  transitions  involving  ground 
state  (or  stretching)  and  /-doubled  bending  vibrational  states  [30], 


Ground,  stretching 

Branching 

1  -doubled 

Branching 

transition 

ratio 

transition 

ratio 

Pk(J ) 

j2  K 2 

rPkU ) 

(K—J)(K—J+l) 

47(27+1) 

7(27+1) 

PPk(J ) 

(-K-J)(-K-J+l) 

47(27+1) 

QkU) 

K2 

rQkU) 

(7-tf)(7+AT+l) 

47(7+1) 

7(7+1) 

pQk(J) 

(J+K)(J—K+1) 

47(7+1) 

rRk{J ) 

(J+K+l)(J+K+2) 

RkU ) 

(X+l)2-*2 

4(7+l)(27+l) 

(7+l)(27+l) 

pRk(J) 

(7— JST+1X7— JST+2) 
4(7+l)(27+l) 

Table  V  lists  only  those  coincidences  in  12CH3F  and 
13CH3F  whose  overlap  with  the  CCU  laser  is  within  one 
atmospheric  pressure-broadened  Lorentzian  linewidth  of  the 
rovibrational  center  frequency.  Of  the  new  rovibrational  coin¬ 
cidences  that  arise,  many  involve  CO2  lines  that  experience 
no  low-pressure  coincidences,  whereas,  the  rest  represent 
additional  coincidences  with  CO2  lines  already  known  to 
have  a  low-pressure  coincidence.  In  both  cases,  atmospheric 
pressure  broadening  of  the  rovibrational  transitions  often 
causes  a  given  laser  line  to  excite  many  rotational  states 
simultaneously.  The  four  rotational  frequencies  reported  in 
Table  V  correspond  to  the  rovibrational  transition  whose 
(Jp,Kp,Vo)  — »■  ( J'p,K'p,Vj )  transition  most  closely  coincides 
with  the  strongest  modulated  terahertz  absorption. 

B.  Differential  absorption 

The  strength  of  the  pump-induced  THz  DR  signatures 
may  be  calculated  from  the  standard  equation  for  molecular 
absorption  [24,27], 

87T3  t 

a  =  — — 1 ){nL  - nu)\ix\-S(v,v0 ) 
ihc 

=  3.24  x  10-6-^—  (fL  —  fu)0  l(I/|/u|£)|2cm_1,  (1) 

Ay 

where  v  is  the  frequency  (GHz)  of  the  transition  and 
rij  =  2.45  x  1013  fi  6  cnD3  is  the  population  of  the 
molecular  states  connected  by  the  probe  and  pump  (i  = 
L,U)  for  a  nondegenerate  fractional  state  population  fi  and 
concentration  9  ppm  at  300  K  in  an  ambient  atmosphere 
of  760  Torr.  The  factor  containing  /z,  is  the  dipole  matrix 
element  for  this  transition  (debye),  and  S^vo)  is  the  line- 
shape  function  that  simplifies  to  1/tt  Ay  (GHz  1 )  at  the  peak 


FIG.  3.  (Color  online)  Energy-level  diagram  depicting  black 
dashed  arrow:  the  IR  transition,  and  red  and  blue  solid  arrows: 
THz  rotational  transitions  connecting  the  pumped  states  JP,  KP, 
and  J'p,K'p.  The  two  pairs  of  THz  DR  signatures  induced  by  the  IR 
pump  are  illustrated  for  the  dotted  red  line:  V,  and  dashed  blue  line: 
V0  vibrational  levels,  the  sum  of  which  for  the  three  types  of  IR 
transition  ( P ,  0,  and  R  types)  produces  the  right:  measurable  DR 
signature. 

of  a  pressure-broadened  Lorentzian  line  shape  [28,29].  The 
absorption  ujh,  at  yTHz  arises  from  rotational  transitions 
mediated  by  the  permanent  dipole  moment  of  the  molecule 
(/z,  =  /z),  whereas  the  absorption  «/«  at  v/s  arises  from 
rovibrational  transitions  mediated  by  the  appropriate  dipole 
moment  derivative  along  the  respective  normal  coordinate  (£?,■) 
for  the  corresponding  excited  vibrational  level  f /x  =  d/i/d Qj) 
[28], 

Note  that  the  absorption  coefficient  is  proportional  to  the 
square  of  the  appropriate  dipole  moment  mediating  the  tran¬ 
sition,  multiplied  by  the  rotational  transition  branching  ratio 
0  <  |([/|/z|L}|2  <  1  given  in  Table  III  [30].  For  IR  transitions, 
the  dipole  derivative  can  be  a  factor  of  10  (stretching  modes) 
to  100  (bending  modes)  times  weaker  than  the  dipole  moment 
itself,  an  effect  that  will  tremendously  affect  DR  spectroscopy. 
Specifically, /z  =  1.86Dandr//x/ii0;  =  0.2756  D  for  the  CH3F 
V3  transition,  whereas, /x  =  1.90D  and  dix/dQe,  =  —  0.0388  D 
for  the  CH3CI  Ve  transition  [29-31].  Consequently,  the  IR 
transition  absorption  strength  is  approximately  2  orders  of 
magnitude  stronger  in  the  stretching  V3  mode  than  in  the 
bending  V(,  mode.  This  difference  arises  because  of  the  type  of 


TABLE  IV.  Optically  pumped  far-infrared  laser  coincidences  for  methyl  fluoride  and  methyl  chloride. 


Molecule 

Laser  line 

IR  transition 

IR  hansition 
offset  (MHz) 

*o" 

R;  (GHz) 

Rf  (GHz) 

Rj  (GHz) 

12CH3F 

9P(20) 

02(12) 

39.78  [25] 

612.409 

663.365 

604.297 

654.582 

13CH,F 

9P(32) 

R3(4) 

-24.25  [25] 

198.855  [36] 

248.559  [36] 

245.350  [36] 

294.406  [36] 

12CH,  35C1 

9P(26) 

*03(6) 

20  [26] 

159.477 

186.050 

158.894 

185.370 

12CH3  35C1 

9R(12) 

rR6(U) 

-30  [26] 

292.190 

318.733 

317.569 

344.009 

12ch3  37C1 

9P(38) 

RP  3(13) 

-50  [26] 

340.093 

366.227 

312.806 

338.850 
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TABLE  V.  Atmospheric  pressure  DR  parameters  for  12CH3F  and  13CH3F. 


Laser  line 

IR  transition 

IR  transition 
offset  (GHz) 

Ro 

vj  (GHz) 

Au(Vj)le  (km^1  kj^1) 
n(Vj) 

< 

vj  (GHz) 

A a(Vj)/e  (km-1  kJ_1) 

n<G) 

K 

Vj  (GHz) 

A a(Vj)ls  (knU1  kJ_1) 

n  {vj) 

Rt 

Vj  (GHz) 

Aa{Vj)ls  (knT1  kj_1) 

H(Vj) 

12CH3F 

9R(18) 

2i(l) 

-2.135  {l}a 

102. 14 1  { 1  }a 

100.786{  l}a 

-4.43 

4.38 

0.27 

0.25 

9R(20) 

2 1:7  (12) 

0.042  {2}a 

612.072  {6}a 

663.000  (6}a 

603.900  {6}a 

654.152  (6)a 

61.03 

-69.14 

-60.22 

68.22 

3.06 

3.04 

3.05 

3.04 

2b(13) 

-1.198  { 1 3 } a 

712.291  { 1 3  }a 

702.391  { 13}a 

-1.61 

1.59 

1.03 

0.80 

Rod) 

1.585  {0}a 

51.072  (0}a 

102.143  (0}a 

50.395  (0)a 

7.64 

-10.18 

7.54 

1.90 

1.00 

1.93 

9P(22) 

2l4:16(18) 

0.0007  (15}a 

914.385  (15}a 

965.016  (15)a 

901.451  (15}a 

951.385  { 15}a 

1.85 

-2.41 

-1.83 

2.37 

1.77 

1.74 

1.74 

1.72 

13CH3F 

9R(08) 

Ro:  16(20) 

0.535  (10}a 

992.507  {3}a 

1041.930  (3}a 

1028.550  {3}a 

1077.320  (3}a 

300.50 

-315.75 

-311.70 

326.72 

4.49 

4.62 

4.60 

4.46 

9P(32) 

Ro:4(4) 

-0.025  (3}a 

198.887  (0)a 

248.598  {0}a 

245.394  {0ja 

294.458  (0}a 

112.30 

-136.47 

-134.71 

158.70 

2.92 

4.40 

4.38 

3.55 

“The  K  quantum  number  of  the  transition  closest  to  the  CO2  laser  line  is  listed  in  curly  brackets  beside  the  IR  transition  offset  frequency.  The 
K  quantum  number  of  the  transition  with  the  lowest  n(v)  is  listed  in  curly  brackets  beside  the  THz  DR  frequencies,  and  the  values  of  Aa/e 
and  n(v)  are  only  listed  for  these  transitions. 


motion  involved,  not  because  of  the  differing  masses  of  the  two 
molecules,  suggesting  a  more  universal  insight:  DR  signatures 
are  generally  stronger  in  stretching  modes  than  bending  modes. 
Unfortunately,  molecules  tend  to  have  more  bending  modes 
than  stretching  modes,  which  can  be  challenging  if  those 
modes  span  the  region  over  which  a  CO2  laser  can  be 
tuned. 

The  strength  of  the  pump-induced  THz  DR  signature  may 
now  be  estimated.  For  an  optically  thin  cloud  (a-\  Hyd  <  1) 
diluted  in  air  to  concentration  6  with  IR  optical  depth  UiRd  < 
1,  the  total  number  of  photons  absorbed  (i.e.,  the  number  of 
molecules  photoexcited  per  pulse)  is 

s 

A  pumped  —  .  OCjRd,  (2) 

hvIR 

where  e  is  the  pump  energy  per  pulse  (J),  d  is  the  path  length 
of  the  laser  in  the  cloud,  and  a/R  can  be  calculated  from  Eq.  ( 1 ) 
with  the  assumption  n (/  ~  0.  From  this,  it  is  easy  to  show  that 
the  fraction  of  the  ul  molecules  pumped  per  laser  pulse  in  a 
beam  of  radius  rjR  (m)  is 


j.  _  ea/ff 

W d~  (itrjR)(hVlR)nL 


6.35 


mm\2z 

A  vIRrjR 


(3) 


where  the  transition  matrix  for  this  rovibrational  transition  is 
the  product  of  the  square  of  the  corresponding  dipole  derivative 
(debye)  and  the  appropriate  transition  branching  ratio  Table  III. 


These  calculations  assume  there  is  no  pump  saturation,  but 
for  highly  intense  pump  beams,  saturation  effects  must  be 
considered. 

If  the  DR  signature  is  observed  on  a  time  scale  shorter  than 
the  collisional  relaxation  time  (~100  ps  at  atmospheric  pres¬ 
sure),  then  the  only  rotational  transitions  with  photoinduced 
population  changes  are  Rq,Rq,R+,  and  /?“  (Fig.  3).  The 
pump-induced  change  in  the  absorption  coefficients  (1)  arises 
from  a  change  in  the  associated  population  differences  from 
equilibrium  ( nR  —  ny)  to  pumped  ( n y  ±  An  —  ny)  values, 
where  An  =  Npumped/W  jR  and  V7 R  —  Tcr2IRd  is  the  volume 
occupied  by  the  IR  beam  inside  the  gas.  The  measured  THz 
DR  spectrum  is  the  pump-induced  change  in  the  rotational 
absorption  coefficient, 


Aa  —  ± 


87T 


2  r 


3h2c 


vTHz 

AVtHz 
Avir 


\(U\n\L)\2 


01 1 R 


Air*  +  (v/r  VpUmp)2  V/r  (tt fjR j 


(4) 


derived  from  Eqs.  (1)  and  (2)  as  a  product  of  terms  that 
depend  on  the  rotational  transition  (first  bracket)  and  on 
the  rovibrational  transition  (second  bracket),  where  vTHz 
corresponds  to  the  monitored  rotational  transition  with  the 
indicated  transition  matrix  element.  Notice  that  a  Lorentzian 
line-shape  factor  has  been  added  to  account  for  the  reduced 
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IR  absorption  caused  by  the  offset  (reported  in  Table  V) 
between  the  pump  and  the  atmospheric  pressure-broadened 
rovibrational  transition  frequencies.  (Figure  2  indicates  the 
pressure  at  which  the  Lorentzian  line  shape  must  be  replaced 
by  a  Gaussian.)  Because  A  a  grows  linearly  with  IR  pump 
energy  (e)  per  pulse,  the  DR  absorption  strengths  are  reported 
as  A  ale  in  units  of  kj  1  km-1  in  Table  V. 


C.  Overlap  enhancement 

Since  the  THz  rotational  transitions  undergo  the  same 
pressure  broadening  as  the  IR  rovibrational  transitions,  the 
THz  DR  spectra  of  nearby  pumped  rotational  transitions  may 
also  overlap.  This  overlap  may  strengthen  or  may  weaken  the 
THz  DR  signal  from  the  low-pressure  value  Eq.  (4)  in  which 
only  one  rovibrational  transition  is  pumped.  The  multiplicative 
factor  by  which  a  given  transition  is  strengthened  or  is 
weakened  is  significantly  less  than  the  number  of  overlapping 
transitions  because  the  constituent  line  strengths  may  have 
differing  degeneracies  and  thermal  populations  as  well  as 
overlaps  with  the  opposite  sign  of  A aj.  (For  a  given  laser 
line,  the  overlapping  transitions  in  CH3F  are  labeled  by  the 
subscripts  in  the  “IR  transition”  column  of  Tables  IV  and  V.) 
The  overlap  enhancement  factor  in  Table  V, 


n (vj)  = 


£" 


(AM2 


KAwp)2 


A  °tj 


(5) 


estimates  how  many  times  stronger  or  weaker  the  DR  signal 
A ctj  will  be  at  frequency  Vj  because  of  the  overlap  from  nearby 
DR  features  of  strength  Aor,  located  at  u, .  Here,  the  summation 
is  over  all  n  DR  signatures  associated  with  all  IR  transitions 
excited  by  a  given  laser  line  where  the  A  a,  and  Aaj  are 
calculated  in  Eq.  (4).  When  the  THz  transitions  overlap,  the 
largest  A aj  necessarily  produces  the  smallest  value  of  n(t>) 
at  the  transition  labeled  {/?},  so  this  is  the  value  provided  in 
Table  V.  The  peak  of  the  summed  DR  signature  will  be  near 
V(r)  ,  so  this  is  the  most  appropriate  frequency  to  use  for  the 
THz  probe. 

The  product  AajTl(Vj)/s  represents  the  final  estimate  of 
the  THz  DR  strength  for  a  given  IR-THz  combination  in  the 
specificity  matrix.  This  product  can  be  used  to  ascertain  which 
of  the  old  or  new  coincidences  are  most  promising  for  remote 
sensing  for  a  given  isotopomer.  Remember  that  AajYl(Vj), 
which  may  be  positive  or  negative,  represents  a  pump  energy- 
dependent  (e)  increase  or  decrease  in  the  ambient  THz  absorp¬ 
tion  strength  as  q'thz('2')  “THz (Vj)  +  AajTl(Vj).  Although 
usually  observed  as  a  short-lived  pump-induced  change  in 
absorption  strength,  THz  emission  can  actually  occur  when 
—  Aar;-n(vy-)  >  «THz(v/).  Because  atmospheric  collisional 
redistribution  of  population  is  minimal  on  time  scales  shorter 
than  100  ps,  AajYl(Vj)/e  only  depends  on  the  absorption 
strength  of  the  IR  rovibrational  molecular  transition,  the 
detuning  of  the  spectral  overlap  with  the  nearest  CCU  laser 
line,  and  the  aggregate  strength  of  the  overlapping  rotational 
transitions  involving  states  simultaneously  excited  by  the 
laser  [n(v;)].  Note  that  the  absorption  strength  of  the  IR 
rovibrational  transition  depends  on  the  type  of  vibrational 
mode  for  a  given  molecule  (dipole  derivative  for  stretching 
or  bending  modes),  the  specific  rovibrational  transition  type 


involved  (branching  ratios  for  P,  Q,  or  R),  and  the  rotational 
partition  function  (or  more  precisely,  the  ground  rotational 
state  population  riL  of  the  pumped  rovibrational  transition). 
Not  surprisingly,  the  optimal  rovibrational  transitions  for  DR 
spectroscopy  are  those  with  the  strongest  IR  absorption,  i.e., 
those  nearest  the  peak  of  the  IR  spectrum. 

Finally,  it  has  already  been  pointed  out  that  the  pressure- 
broadened  Lorentzian  line-shape  function  falls  off  less 
abruptly  than  the  Doppler-broadened  Gaussian  line-shape 
function,  so  overlaps  are  more  readily  excited.  The  decision  to 
include  only  those  overlaps  within  one  line  width  A  Vjp  of  the 
rovibrational  line  center  was  arbitrary.  Recall  the  calculation 
assumed  the  CO2  laser  was  operating  continuously;  however, 
if  the  laser  was  pulsed,  the  increased  spectral  bandwidth  of 
the  laser  pulses  should  also  be  considered  in  the  coincidence 
calculation.  For  the  remote  sensing  application,  pump  pulses 
of  durations  comparable  to  the  hard-shell  gas  kinetic  collision 
time  (~100  ps)  are  required,  corresponding  to  a  spectral 
bandwidth  of  10  GHz  that  is  (not  coincidentally)  comparable  to 
the  pressure-broadened  linewidth  Avp  of  the  IR  rovibrational 
transition.  Clearly,  even  more  coincidences  may  be  excited 
when  such  short  pulses  are  used,  but  the  pump  efficiency  is 
reduced  by  the  fraction  of  pump  radiation  that  falls  outside  the 
spectral  bandwidth  of  the  IR  absorption  line. 


III.  SPECIFIC  EXAMPLES 

Having  outlined  the  basic  principles  underlying  the  exci¬ 
tation  of  a  DR  signal,  we  next  explore  the  variety  of  DR 
signal  types,  their  strengths,  and  how  they  depend  on  the 
type  of  photoexcited  rovibrational  transition  (Fig.  3).  Again, 
we  restrict  ourselves  to  photoexcitation  on  time  scales  short 
compared  to  the  fastest  collisional  relaxation  time  of  the 
molecule  (~  1 00  ps),  so  only  rotational  transitions  R([ ,  RQ  ,  Rf , 
and  R-  are  involved. 


A.  Methyl  fluoride 

We  start  with  the  simpler  case,  CH3F,  for  which  the 
stretching  mode  V3  of  both  12CH3F  and  13CH3F  isotopomers 
is  photoexcited  by  a  CO2  laser.  Considering  first  12CH3F, 
the  9  R  (20)  line  of  a  CO2  laser  is  known  to  photoexcite  the 
<22(12)  [and  perhaps  the  Q i(12)]  transition  at  low  pressure. 
The  resulting  DR  signature  has  components  at  612.409 
(Ro  ),  663.365  (R+),  604.297  (Rj),  and  654.582  GHz  (R+). 
The  absorption  decreases  (A aj  <  0)  in  transitions  R( j"  and 
R 3  and  increases  (A ctj  >  0)  in  transitions  R^  and  /£ 
[Figs.  3  and  4(a)],  illustrating  that  a  characteristic  Q- branch 
THz  DR  signature  consists  of  four  spectral  features  in  the  form 
of  conjugate  doublets  (two  R  and  two  R+)  separated  by  ~2 B 

5 1  GHz.  Assuming  there  is  no  IR  pump  saturation,  Eq.  (4) 
indicates  the  strengths  of  these  four  DR  features  are  linearly 
proportional  to  the  pump  pulse  intensity,  their  THz  frequency, 
and  the  associated  branching  ratios  (Table  III).  Each  doublet 
consists  of  the  identical  rotational  transitions  in  Vo  and  V3 
separated  in  frequency  by  ~2(Z?3  —  Bq)  U  +  1),  a  separation 
of  8.1  GHz  for  J  =  11  -»  12  and  8.8  GHz  for  J  =  12  ^ 
13.  Given  that  Avp  —  2.3  GHz,  these  four  features  are  well 
resolved  at  atmospheric  pressure. 
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Frequency  (GHz) 


FIG.  4.  (Color  online)  (a)  Aggregate  DR  signature  for  the  I2CH3F 
transitions  pumped  by  the  9/’(20)  laser  line,  plotted  above  the 
constituent  labeled  components  from  Table  V  (to  scale).  Energy-level 
diagram  depicting  both  infrared  pump  transitions  and  rotational 
transitions  connected  to  the  pumped  states,  (b)  Aggregate  DR 
signature  for  the  CH3  35  Cl  transitions  pumped  by  the  9P  (26)  laser 
line,  plotted  above  the  constituent  components  (to  scale)  labeled 
center.  Energy-level  diagram  depicting  both  infrared  pump  transitions 
and  rotational  transitions  connected  to  the  pumped  states.  The 
transitions  PR0( 24),  RR0(24),  and  RPj( 31)  from  Table  VI  are  not 
depicted  because  they  occur  at  higher  frequencies. 

As  noted  above,  the  atmospheric  pressure -broadened  width 
of  the  rovibrational  transitions  allows  a  given  laser  line  to  pho- 
toexcite  many  neighboring  transitions,  primarily  those  transi¬ 
tions  with  differing  K  values,  further  increasing  the  strength 
of  the  DR  signature.  In  12CH3F,  the  9P(20)  laser  line  now  co¬ 
incides  with  nine  additional  atmospheric  pressure-broadened 
rovibrational  transitions  [0i(12)  —  07(12),  0  i3(  1 3),  and 
Pq(  1 )] .  When  a  series  of  transitions  with  constant  J  but  varying 
K  are  pumped  [e.g.,  0 i (12)  —  07(12)],  it  is  listed  as  a  single 
entry  in  Table  V  with  {R}  =  {K }  for  the  transition  with  smallest 
n(Vj).  Figure  2  shows  how  the  multiplier  V\(vK )  for  each 
of  the  Qk(12)  transitions  ( K  —  1-7)  grows  with  increasing 
pressure,  and  Figs.  4(a)  and  5(a)  show  how  their  aggregate 
DR  signature  is  the  sum  of  all  the  constituent  overlapping 


Frequency  (GHz) 

FIG.  5.  Calculated  DR  signatures  for  all  laser  concidences  with 
(a)  12CH3F  and  (b)  13CH3F.  The  coincident  CCF  laser  lines  are  listed 
on  the  right. 


0-branch  signatures  for  a  portion  of  the  THz  region.  Notice 
that  n(uy)  =  1  for  02(12)  until  the  pressure  broadening 
exceeds  the  Doppler  broadening  of  the  IR  transition,  after 
which  there  is  a  rapid,  often  nonmonotonic  growth  in  II  (y,-) 
for  all  transitions  as  pressure  increases.  Notice  Il(y,)  for 
{7C}  =  6  only  reaches  3.04  at  atmospheric  pressure,  not 
the  theoretical  maximum  of  6.15  if  all  seven  K  levels  were 
equally  populated  and  were  equally  pumped.  This  reduction 
occurs  because  of  two  effects:  the  unfavorable  branching  ratios 
for  low -K  0-branch  transitions  and  the  greater  frequency 
separation  of  the  rovibrational  transitions  in  the  IR  than  the 
associated  rotational  transitions  in  the  THz.  For  these  same 
reasons,  all  other  transitions  have  larger  Il(v/)  values  because 
their  isolated  DR  signature  Aoq  is  weaker.  In  other  words,  at 
atmospheric  pressure,  the  rovibrational  transition  that  makes 
the  largest  contribution  to  A  a  is  06(12),  not  02(12)  as  at  low 
pressure  [Fig.  4(a)]. 

At  atmospheric  pressure,  the  9P(20)  line  also  excites  other 
rovibrational  transitions  besides  the  Qk(  12)  series,  namely, 
013(13)  and  Po(l)  [Fig.  5(a)],  These  are  quite  weak  because  of 
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9P(xx)  1  OP(xx) 


Frequency  (GHz)  Frequency  (GHz) 

FIG.  6.  Calculated  DR  signatures  for  all  laser  concidences  with  CH3  35C1.  The  coincident  CO2  laser  lines  are  listed  on  the  right. 


the  low  thermal  population  in  the  respective  (Jp,Kp,V 0)  states. 
Laser  lines  adjacent  to  9P(20)  also  produce  new  coincidences, 
but  the  DR  signatures  are  weaker  [Fig.  5(a)].  This  suggests  a 
trend  that  becomes  more  obvious  when  CH3CI  is  studied:  The 
overlap  of  rovibrational  transitions  and  CCL  laser  lines  can  be 
likened  to  a  moire  pattern  formed  by  two  overlapping  grids 
with  slightly  varying  spatial  frequencies.  Overlaps  more  often 
come  in  clusters  of  CO2  laser  lines  instead  of  through  random 
isolated  lines. 

P-  and  /(-type  rovibrational  transitions  also  manifest  four 
THz  DR  components,  but  these  occur  as  a  triplet  composed 
of  a  central  doublet  flanked  by  singlets  of  the  opposite  sign 
(Fig.  3).  The  differing  shapes  of  the  P-,  Q-,  and  /(-dependent 
DR  spectra,  particularly,  the  associated  A oq  sign  differences, 
represent  an  additional  degree  of  recognition  specificity  at  a 
specific  THz  probe  frequency.  A  typical  example  of  a  R -branch 
DR  signature  comes  from  the  A’dT)  transition  in  13CH3F, 
optically  pumped  by  the  9P(32)  line  of  a  CO2  laser  [Fig.  5(b)]. 
A  P -branch  DR  signature  would  be  identical  except  flipped  in 


sign.  The  central  doublet  consists  of  two  components  of  the 
same  sign  separated  by  v  ~  2(Z?3  —  Bo)  (J  +  1)  =  3.4  GHz, 
whereas,  the  singlets  of  opposite  sign  are  separated  from  the 
doublets  by  v  ~  2B  &  50  GHz.  Consequently,  the  DR  signature 
is  spread  twice  as  far  in  frequency  as  0-branch  DR  signatures. 
The  central  doublet  is  barely  resolved  and  very  strong  at 
atmospheric  pressure. 

It  is  striking  how  sparse  the  THz  DR  spectra  are  in 
Figs.  5(a)  and  5(b).  Although  12CH3F  and  13CH3F  have 
13  and  22  pressure-broadened  coincidences,  respectively,  most 
of  them  occur  in  clusters  of  adjacent  rovibrational  transitions 
of  constant  J  but  varying  K  excited  by  the  same  laser  line. 
Consequently,  there  are  only  five  and  two  distinguishable  coin¬ 
cidences  for  12CH3F  and  13CH3F,  respectively,  each  of  which 
has  up  to  four  frequency-resolvable  THz  DR  components.  No 
CO2  laser  line  that  excites  one  isotopomer  excites  the  other. 
There  are  no  DR  signatures  of  one  isotopomer  within  the 
same  2  A  Vp  ^  5 -GHz  window  of  any  DR  signature  from  the 
other  isotopomer.  Indeed,  the  recognition  specificity  matrices 
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FIG.  7.  Calculated  DR  signatures  for  all  laser  concidences  with  CH3  37C1.  The  coincident  CO2  laser  lines  are  listed  on  the  right. 


of  12CH3F  and  13CH3F  are  so  sparse  that  the  isotopomers  are 
easily  distinguished  through  the  selection  of  the  CO2  laser  line 
and  corresponding  THz  probe  frequency. 

B.  Generalizations  to  other  molecules 

To  explore  how  general  a  result  this  is  for  stretching 
mode  DR  spectra,  consider  more  massive  or  complex  pro¬ 
late  molecules  and  their  correspondingly  smaller  rotational 
constants.  Because  individual  absorption  strengths  a-I  H/  will 
weaken  with  increasing  molecular  mass  and  rotational  parti¬ 
tion  function,  the  correspondingly  reduced  tii  and  am  might 
also  be  expected  to  weaken  the  DR  signals  Aa;-  n(v;-)/e. 
However,  the  number  of  rovibrational  IR  coincidences  and 
rotational  state  overlaps  [TTC )]  increase  with  increasing 
molecular  mass  and  complexity  because  the  rotational  spac- 
ings  shrink  and  the  peak  THz  absorption  redshifts.  These 
effects  partially  or  completely  compensate  for  the  increased 
partition  function  and  reduced  n  / .  The  increasing  number  of 


coincidences  with  CCU  laser  lines  will  also  add  more  elements 
to  the  recognition  specificity  matrix,  permitting  selection  of 
favorable  DR  features  in  atmospheric  transmission  windows. 
Nevertheless,  the  specificity  matrix  is  very  large,  containing 
approximately  60  x  200  x  2  =  24  000  unique  combinations 
of  sixty  CO2  laser  frequencies,  two  hundred  5-GHz-wide  win¬ 
dows  from  0  to  1  THz,  and  ±  sign  of  A  a  j ,  respectively.  Only 
in  the  rarest  of  cases  when  two  molecules  excited  by  the  same 
laser  line  and  produce  a  DR  signal  of  the  same  sign  in  the  same 
5-GHz-wide  window  will  this  affect  the  distinguishability  of 
molecules  through  their  DR  signatures.  Such  false  positive 
detections  can  be  avoided  by  simply  photoexciting  a  different 
coincidence. 

The  challenge  of  multiple  coincidences  requires  a  strategy 
for  identifying  the  strongest  DR  signature.  Several  insights 
from  the  analysis  of  CH3F  arise.  All  isotopomers  (and  their 
isotopic  abundances)  must  be  considered  when  identifying  the 
strongest  DR  signature  for  a  given  molecule.  In  the  case  of 
methyl  fluoride,  the  largest  AajYl(vj)/s  for  the  photoexcited 
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/(-branch  rovibrational  transitions  of  13CH3F  is  seven  times 
larger  than  the  strongest  0 -branch-excited  values  in  12CH3F. 
Only  13C’s  rarity  (1%  isotopic  abundance)  prevents  13CH3F 
from  being  the  preferred  isotopomer. 

Next,  consider  the  factor  \{U\j^\L)\2  in  AajTl(Vj)/e  for 
these  rovibrational  transitions,  which  may  be  expressed  as  the 
product  of  (d/u-ldQ)2  and  the  quantum-mechanical  branching 
ratios  in  Table  II.  Because  of  the  latter,  (7-branch  transitions 
are  only  favored  when  /  &  K.  However,  levels  with  ./  ~  K 
have  the  weakest  thermal  rotational  population,  and  these  same 
branching  ratios  dramatically  weaken  the  corresponding  AJ  = 
1  DR  signatures.  For  a  given  K ,  the  branching  ratio  for  0- 
branch  transitions  tends  to  zero  as  J  increases,  whereas,  P- 
and  A1 -branch  transitions  tend  to  Because  low- AT  transitions 
have  closer  THz  spectral  spacing  and  comparatively  larger 
thermal  populations  than  high -K  transitions,  it  is  evident  that 
photoexciting  low-Al  P-  or  /(-branch  transitions  will  yield  the 
largest  n«  enhancement  factors. 

Other  factors  also  favor  P-  and  R  -branch  transitions 
over  the  0  branch.  It  has  already  been  mentioned  that 
the  strongest  DR  signatures  arise  when  transitions  near  the 
peaks  of  the  IR  spectra  are  excited.  For  increasingly  massive 
molecules,  the  rotational  quantum  numbers  J  associated  with 
these  peaks  will  increase,  and  the  rotational  constants  B 
and  their  differences  ( B,  —  Bo)  will  decrease.  Now,  consider 
again  the  2(B^  —  Bq)(J  +1)  spacings  of  the  doublets  in  the 
respective  transitions  (ignoring  higher-order  distortion  terms) 
as  illustrated  in  Fig.  3.  In  general,  features  are  resolvable  when 
2(S3  —  B0)(J  +  1)>  A  Vp.  In  P- and  R  -branch  transitions,  the 
doublet  components  have  the  same  sign  of  Aa,  so  decreasing 
(Z?3  —  Bo)  or  J  simply  increases  the  overlap  and  effectively 
doubles  the  strength  of  the  central  doublet.  By  contrast,  the 
doublet  components  in  0  -branch  transitions  have  opposite 
signs,  so  decreasing  (B3  —  Bo)  or  J  can  weaken  the  doublet 
to  the  point  of  disappearance.  This  explains  why  (7-branch 
transitions  often  have  DR  signatures  for  which  FI(v;)  is  less 
than  1  and,  more  importantly,  suggests  that  P-  and  R -branch 
transitions  are  more  resilient  than  (7-branch  transitions  to  the 
trade-off  between  decreasing  (B,  —  Bo)  and  increasing  J. 

Summarizing  our  findings  to  this  point,  the  strongest 
DR  signatures  involving  vibrational  stretching  modes  will 
be  obtained  when  photoexciting  rotational  transitions  with 
low  K  whose  J  values  correspond  to  P-  or  /(-branch 
rovibrational  transitions  near  the  peak  of  the  corresponding 
IR  spectra.  Although  these  rules  were  derived  for  symmetric 
top  molecules,  their  extension  to  the  stretching  modes  of  other 
types  of  molecules  is  straightforward  and  is  simplified  by  the 
fact  that  many  of  the  asymmetric  rotor  distortions  and  hyper- 
fine  splittings  are  unresolved  within  the  pressure-broadened 
linewidth  of  the  DR  signal.  These  transitions  constitute  the 
“overlap  enhancement  factor”  n(v,)  for  asymmetric  rotors. 

C.  Methyl  chloride 

Turning  now  to  CH3CI,  the  calculated  THz  DR  spectra 
in  Figs.  6  and  7  reveal  that  both  isotopomers,  CH3  35C1  and 
CH3  37 Cl,  exhibit  many  more  coincidences  than  in  CH3F. 
Although  it  is  true  that  the  B  rotational  constant  is  roughly  half 
that  in  CH3F  (13.3  GHz  for  CH3C1,  25.5  GHz  for  CH3F),  that 
fact  alone  should  only  double  the  number  of  DR  signatures 


over  CH3F.  Instead,  CH3  35C1  and  CH3  37C1  have  154  and 
140  coincidences  within  one  atmospheric  pressure-broadened 
linewidth,  respectively,  seven  to  ten  times  more  than  CH3F. 

The  difference  is  caused  by  the  fact  that  the  laser  is  now 
coincident  with  the  14  vibrational  bending  level  rather  than  the 
V3  stretching  level.  Because  bending  modes  (/  =  ±1)  couple 
vibrational  and  rotational  (K)  motion  to  produce  /-doubled 
spectra,  the  selection  rules  and  branching  ratios  in  Table  III 
produce  much  richer  rovibrational,  rotational,  and  THz  DR 
spectra  (Figs.  1,  6,  7).  In  particular,  the  ±K  degeneracy  is 
removed  by  large  energies  that  are  proportional  to  Kl,  and 
rovibrational  transitions  connect  states  of  quantum  number 
K  in  Vo  with  states  of  quantum  number  (K  —  l)  in  14- 
Besides  doubling  the  number  of  rotational  and  rovibrational 
transitions,  /  doubling  introduces  large  energy  and  frequency 
splittings  in  14,  so  fewer  rovibrational  transitions  occur 
within  an  atmospheric  pressure-broadened  linewidth.  Indeed, 
atmospheric  pressure  broadening  never  adds  coincidences  in 
CH3CI  with  adjacent  rovibrational  transitions  that  preserve  J 
but  change  K  for  a  given  laser  line,  in  stark  contrast  with  CH3F. 

Instead,  it  is  often  found  that  a  given  laser  line  photoexcites 
a  complex  combination  of  overlapping  and  nonoverlapping 
P-,  0-,  and/or  /(-branch  transitions  simultaneously,  typically 
involving  widely  separated  rotational  states  and  THz  transition 
frequencies  (see  Figs.  6  and  7).  The  DR  signature  can  be  under¬ 
stood  as  a  collection  of  these  individual  P-,  Q-,  and  /(-branch 
DR  signatures  (Fig.  3)  whose  strengths  and  frequencies  vary 
widely  based  on  the  idiosyncrasies  of  the  constituent  rotational 
quantum  numbers,  energies,  and  transition  branching  ratios. 
The  aggregate  DR  spectrum  for  a  given  laser  line  is  quite 
unique,  providing  additional  specificity. 

For  example,  consider  the  9/426)  laser  line,  which  has  a 
known  low-pressure  coincidence  with  the  *03(6)  rovibrational 
transition  in  CH3  35C1  [Fig.  4(b)].  At  atmospheric  pressure, 
nine  more  P-,  Q-,  and  /(-branch  rovibrational  transitions  are 
photoexcited,  including  the  entire  0  -branch  sequence  from 
SG3(4)  to  s03(8),  plus  */(2( 7),  *P4(8),  rR0(24),  pR0(24),  and 
K Pj(3 1 ),  producing  DR  spectra  that  span  106-848  GHz.  Inter¬ 
estingly,  the  strongest  of  these  features  are  not  those  associated 
with  the  known  low  pressure  *03(6)  coincidence  but  with  the 
new  */((|(24)  coincidence.  A  contributing  reason  for  this  is 
the  increased  overlap  of  the  opposite  sign  0-branch  doublets 
spaced  2(Bj  —  Bq)(J  +  1)  apart:  ( Bf]  —  Bo)  rs  49  MHz  in 
CH3CI  while  («3  -  Bo)  s*  339  MHz  in  CH3F.  Although  the 
coincidences  in  CH3CI  typically  involve  states  with  larger 
J  than  CH3F,  atmospheric  pressure  broadening  causes  the 
0-branch  doublets  to  overlap  strongly  enough  to  reduce  n(vy) 
below  1  in  most  cases. 

Comparing  the  relative  strengths  of  the  constituent  transi¬ 
tions  excited  by  the  9/426)  line  provides  helpful  insight  into 
the  subtle  interplay  of  factors  that  contribute  to  the  complex 
DR  signature.  First  of  all,  because  the  laser  frequency  is  within 
one  pressure-broadened  linewidth  of  all  ten  coincident  IR 
transitions,  the  strength  of  the  corresponding  DR  signatures 
depends  primarily  on  the  branching  ratio  of  the  transition. 
Table  III  indicates  that,  as  K  — »•  J,  the  IR  transitions  rR 
and  pP  grow  much  stronger  than  the  p R  and  RP  transitions. 
For  example,  the  branching  ratios  of  the  IR  transitions  R Rid) 
and  *74(8)  are  0.229  and  0.022,  respectively,  explaining  why 
their  contributions  to  the  DR  signature  differ  by  an  order  of 
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TABLE  VI.  Abbreviated  table  of  atmospheric  pressure  DR  parameters  for  CH3  35C1.  The  complete  table  may  be  found  in  Ref.  [39], 


Laser  line  IR  transition 

IR  transition 
offset  (GHz) 

R0 

vj  (GHz) 

Aa(Vj)/e  (km-1  kj_1) 

H(Vj) 

< 

vj  (GHz) 

A a(Vj)/e  (km-1  kj_1) 

n  (vj) 

K 

Vj  (GHz) 

A a(Vj)/s  (km~*  kj_1) 

n  (vj) 

Kt 

Vj  (GHz) 

Aa(i'j)fe  (km-1  kj-1) 

H(Vj) 

9/*  (26)  *03(4) 

1.09 

106.324 

132.902 

132.416 

0.072 

-0.128 

0.072 

1.008 

-0.500 

0.763 

*03(5) 

0.60 

132.902 

159.477 

132.416 

158.894 

0.263 

-0.364 

-0.148 

0.268 

0.243 

-0.150 

-0.371 

0.136 

*03(6) 

0.02 

159.477 

186.050 

158.894 

185.370 

0.533 

-0.669 

-0.394 

0.550 

0.102 

-0.676 

-0.092 

0.724 

*/(2(7) 

-1.28 

186.064 

212.636 

211.858 

238.329 

0.457 

-0.528 

-0.482 

0.557 

0.990 

1.709 

1.825 

0.951 

*03(7) 

-0.66 

186.050 

212.620 

185.370 

211.843 

0.760 

-0.905 

-0.624 

0.787 

0.595 

0.998 

-0.637 

-1.115 

*03(8) 

-1.44 

212.620 

239.187 

211.843 

238.312 

0.815 

-0.942 

-0.709 

0.847 

-1.108 

-0.244 

1.238 

0.630 

*/L(  8) 

0.69 

212.598 

239.162 

185.351 

211.821 

0.041 

-0.049 

-0.024 

0.033 

-22.055 

-4.862 

-16.765 

-26.324 

*Ro(24) 

1.04 

637.057 

663.512 

661.071 

687.418 

2.316 

-2.410 

-2.397 

2.491 

0.970 

1.425 

1.430 

0.971 

*Ro(24) 

1.13 

637.057 

663.512 

661.075 

687.422 

2.244 

-2.335 

-2.323 

2.414 

1.001 

1.471 

1.476 

1.002 

*P7(31) 

-0.31 

821.399 

847.750 

792.123 

818.390 

0.162 

-0.168 

-0.153 

0.159 

1.345 

0.988 

0.986 

1.359 

magnitude  [Fig.  4(b)].  Likewise,  the  rQk(J)  branching  ratios 
(Table  III)  tend  to  zero  as  /  -*  K,  so  inspection  of  the  0- 
branch  series  *03(4)  to  *03(8)  in  Fig.  4(b)  reveals  the  Vg  con¬ 
tributions  to  the  0-branch  doublets  disappear  as  J  decreases. 

Figures  6  and  7  suggest  the  strongest  DR  feature  for 
almost  every  coincidence  derives  from  photoexcitation  of  a  P- 
or  /(-branch  transition.  Although  laser  coincidences  usually 
occur  with  isolated  P-  or  /(-branch  transitions,  sometimes 
a  sequence  of  laser  lines  coincides  with  a  sequence  of  IR 
transitions.  A  prominent  example  involves  the  series  of  laser 
lines  spanning  10F’(26)-10F’(38)  where  coincidences  that 
preserve  K  —  9  for  odd  ./  span  pPg(23)-p  Pg(35),  one  pet- 
laser  line.  Many  shorter  sequences  with  this  same  moire  like 
overlap  pattern  of  preserving  K  but  changing  J  were  found  in 
both  isotopomers. 

By  contrast,  0-branch  transitions  cluster  strongly,  so 
multiple  IR  transitions  may  be  excited  by  a  single  laser  line. 
Unlike  V3  of  CH3F  for  which  the  P-,  Q-,  and  /(-branch 
IR  transitions  are  arranged  in  clusters  that  preserve  J  for 
all  K  [Fig.  1(b)],  /  doubling  in  Vg  of  CH3CI  spreads  out 
the  0-branch  transitions  into  clusters  that  preserve  K  for 
all  J  [Fig.  1(c)].  Consequently,  a  given  CO2  laser  line  may 


excite  many  such  /-changing  /(/-preserving  IR  transitions, 
explaining  how  the  9P(26)  line  excites  *03(4)-*03(8).  When 
this  occurs,  an  additional  reduction  in  the  DR  signature  occurs. 
Recall  that  the  spacing  of  the  transitions  and  the  separation  of 
the  0-branch  components  of  opposite  sign  are  both  ~2B.  As  a 
result,  the  positive  component  of  one  DR  feature  overlaps  the 
negative  component  of  the  preceding  DR  feature  [Fig.  4(b)], 
and  the  summed  0-branch  DR  signature  is  much  weaker  than 
its  constituent  components.  Therefore,  0-branch  transitions 
generally  produce  unfavorably  weak  DR  features  for  bending 
vibrational  modes,  just  as  they  did  in  stretching  vibrational 
modes,  but  for  different  reasons. 

So  atmospheric  pressure  broadening  typically  produces 
overlapping  sequences  of  0-branch  DR  features  with  FI(v;) 

<  1  and  isolated  P-  and  /(-branch  DR  features  with  1  <  n  (vj) 

<  2.  However,  inspection  of  Table  VI  indicates  that  n(vy)  can 
sometimes  be  much  larger.  For  example,  n(vy)  ~  25  for  the 
two  DR  signatures  near  212  GHz  associated  with  the  A'  /J4(  8 ) 
transition.  Such  IT  ( Vj )  >  2  cases  are  not  caused  by  the  summa¬ 
tion  of  overlapping  DR  signatures  excited  by  a  single  IR  tran¬ 
sition  as  is  the  case  for  CH3F.  Instead,  these  are  rare  instances 
in  which  very  different  rovibrational  transitions  excited  by  the 
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same  laser  line  produce  THz  DR  signatures  that  coincidentally 
overlap.  In  the  example  at  212  GHz,  the  weak  ^4(8)  DR 
signatures  overlap  with  DR  signatures  17  and  5  times  stronger 
produced  by  the  S/G(7)  and  RQi(S)  transitions,  respectively. 

The  increased  mass  of  CH3  37C1  lowers  the  vibrational 
energy  of  V, 6,  so  the  DR  spectra  of  CH3  37C1  are  similar  in 
form  but  different  in  detail  to  those  of  CH3 35 Cl.  The  isotopic 
abundances  (76%  and  24%  for  35C1  and  37C1,  respectively) 
suggest  that  CH3  35C1  will  be  the  preferred  isotopomer  for 
remote  sensing.  However,  if  a  DR  signature  from  CH3  37C1 
is  three  times  stronger  than  any  from  CH3 35  Cl,  the  less 
abundant  isotope  could  produce  the  strongest  DR  signature. 
Indeed,  if  we  constrain  ourselves  to  look  at  the  region  of 
greatest  practical  interest  below  300  GHz  and  to  account  for 
isotopic  abundances,  we  find  that  the  strongest  DR  signature 
for  CH3CI  is  the  rR3(  10)  transition  of  CH3  37C1  at  287  GHz 
excited  by  the  9 P(  1 6)  laser  line. 

IV.  CONCLUSIONS 

Our  analysis  has  shown  the  potential  of  IR-THz  DR 
spectroscopy  for  the  remote  recognition  of  a  trace  gas  in 
the  atmosphere  given  sufficient  pump  intensity  and  gas 
concentration.  Discrimination  comes  from  the  two  dimensions 
of  a  sparse  recognition  specificity  matrix  derived  from  (1)  the 
few  coincidences  that  occur  between  CO2  pump  laser  lines 
and  rovibrational  transitions  of  the  analyte  and  (2)  the  unique 
changes  in  THz  absorption  the  pump  coincidences  induce.  The 
strength  of  the  DR  signature  depends  sensitively  on  the  type  of 
rovibrational  transition  excited  by  the  laser  coincidence  and  the 
quantum  and  statistical  mechanics  of  the  involved  transitions. 


Atmospheric  pressure  broadening,  which  smears  out  transi¬ 
tions  detected  with  traditional  THz  spectroscopic  techniques, 
actually  assists  the  IR-THz  DR  technique  by  providing  more 
pump  coincidences  and  more  overlapping  transitions  for  a 
given  coincidence.  These  advantages  allow  an  optimal  DR 
signature  to  be  selected  for  a  given  analyte,  and  it  was  shown 
how  to  calculate  this  from  basic  molecular  parameters. 

Generally  speaking,  the  DR  signatures  for  CH3CI  are  much 
weaker  than  those  for  CH3F:  The  transition  dipole  moments 
are  smaller,  and  the  rotational  partition  function  is  larger, 
so  the  rovibrational  pump  absorption  coefficient  is  weaker. 
Weak  IR  absorption,  particularly  for  molecules  whose  pump 
coincidences  involve  bending  modes,  presents  the  greatest 
challenge  to  DR  spectroscopy,  just  as  it  has  limited  the 
performance  of  many  optically  pumped  far-infrared  laser 
gases.  For  heavier  molecules,  this  problem  may  be  alleviated 
as  the  transition  dipole  moments  tend  to  increase  (e.g., 
dn/dQs  =  0.054  D  for  CH3BQ,  and  increasing  complex 
vibrational  energy-level  structures  allow  some  selection 
in  which  vibrational  mode  the  laser  excites  (e.g.,  a  CO2 
laser  may  excite  the  V4,  V7,  or  3Vg  vibrational  modes  of 
CH3CN)  [32,33]. 
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